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ABSTRACT 

The  kinematic  motion  of  (N  +  1)  geostrophlc  (Bessel) 
vortices  Is  studied.   Initially  the  vortices  are  positioned 
near  a  uniformly  rotating  equilibrium  configuration  which 
consists  of  N  vortices  of  equal  strength  7,  equally  spaced 
on  a  circle  of  radius  a,  and  one  vortex  of  strength  7  at 
the  center  of  the  circle.   One  length  scale  K~  is 
associated  with  each  of  the  vortices,  the  limiting  case 
K  =  0  being  the  classical  logarithmic  vortex.   The  exponential 
stability  of  this  equilibrium  configuration  is  studied  over 
the  range  of  the  three  parameters:   N_>2,  0  <  (Ka)  <  00 
and  -  CO  <  {7^/7)    ^  00  •   The  need  for  a  supplementary 
nonlinear  study  is  discussed. 


il 


1.   Continuum  and  Discrete  Description  of  Geostrophic  Motion. 

The  first  systematic  study  of  geostrophic  vortex 
motion  was  begun  by  H.  J.  Stewart  [1,2]*,  using  primitive 
computers  by  present-day  standards.   By  evoking  some  simple, 
but  rather  convincing,  physical  arguments  he  represented 
the  large-scale  closed  isobaric  systems  of  the  atmosphere  by 
discrete  rectilinear  geostrophic  (Bessel)  vortices.   Our 
paper  is  an  extension  of  this  earlier  work. 

In  more  recent  years  with  the  advent  of  large  electronic 
computers,  formally  more  systematic  derivations  of  the 
approximate  non-linear  equations  describing  geostrophic 
motion  have  been  found,   e.g.  by  Obukhov  [3]  and  Morikawa 
[4].   Meteorologists  have  applied  such  equations  to  short 
range  weather  forecasting  studies  using  continuum  (or 
smooth)  initial  data. 

During  the  past  few  years  we  have  been  studying,  mainly 
by  numerical  methods,  a  ni;miber  of  different  types  of  initial 
value  problems,  i.e.  1)  purely  continuum  data,  2)  a  combina- 
tion of  smooth  and  discrete  (vortices)  data,  and  J>)    purely 
discrete  data.  A   primary  interest  is  to  understand  the 


*The  existence  of  Stewart's  early  work  which  apparently  has 
been  ignored  in  the  meteorological  literature,  was  made 
known  to  one  of  us  by  S.  C.  Lowell. 


non-linear  Interacting  development  of  the  motion.   We  have 
done  more  work  on  1)  and  3)  than  on  2)  which  requires  the 
most  computer  space*.   In  this  paper  we  describe  In  part 
some  of  our  work  on  5) :   specifically,  we  study  the 
Interacting  motion  of  N  rectilinear  geostrophic  vortices  of 
equal  strength  7  initially  distributed  approximately  with 
equal  spacing  near  (or  on)  a  circle  of  radius  a  and  of  a 
single  vortex  of  strength  7  placed  initially  near  (or  at) 
the  center  of  the  circle.   All  N  +  1  vortices  have  the  same 
horizontal  length  scale  K~  where  the  limiting  case  K  =  0 
corresponds  to  potential  (logarithmic)  vortices.   This 
paper  consists  of  two  principal  parts:   I)  linearized 
stability  analysis,  an  eigenvalue  problem,  and  II)  numerical 
computation  of  the  non-linear  motions  in  sufficient  detail 
to  interpret  the  stability  results.   The  study  of  linear- 
ized stability  for  the  special  case  of  K  =  0,  7   -  0  goes 
back  many  years  to  an  early  essay  of  J.  J.  Thomson  [8] 
who  was  Interested  in  developing  a  vortex  model  of  an 
atom  (unsuccessful). 

A   single  rectilinear  geostrophic  vortex  in  an 


*Some  of  that  part  of  2)  related  to  describing  hurricane 
trajectories,  using  a  single  vortex  to  represent  the 
hurricane,  has  been  reported  in  [5]  at  a  Symposium  in 
Tokyo  in  I96O;  this  numerical  work  was  carried  out  with 
E.  Isaacson  and  D.  A.  Levine  [6].   Some  results,  comparing 
the  difference  in  numerical  solutions  between  the  use  of 
smooth  data  and  a  combination  of  smooth  data  interacting 
with  a  single  vortex,  were  presented  at  a  Symposium  in 
Mexico  City  in  I965  [7I. 


imconflned   space    is   given  by   the    stream   function 

(1.1)  ^  =  JLK(»c|r-r|) 

where  the  stream  function  ^  for  any  geostrophlc  motion 
(in  our  sense)  is  related  to  the  velocity  components  by 

(1.2)  ,„.,  =,.|£,  ||, 

which  states  that  this  approximate  two-dimensional  descrip- 
tion of  large-scale  nearly  horizontal  atmospheric  motion  is 
divergence-free;  but  in  addition  (1.2)  states  that  the  flow 
is  geostrophic  (and  hydrostatic,  as  well)  to  this 
order  of  approximation,  i.e.  the  inertia  term  in  the 
horizontal  momentum  equations  is  small  compared  to  the 
pressure  gradient  and  coriolis  terms.   The  distance 
|r  -  r  I  =  [(x  -  X  )   +  (y  -  y  )  ]    where  the  vortex 
position  (x  ,y  )  may  be  a  fimction  of  time  if  there  are 

other  elements  in  the  field  of  motion.   K  is  the  Bessel 

o 

function  of  the  second  kind  (zeroth  order),  singular  at 
the  origin. 

For  the  problem  which  we  consider  with  (N  +  1)  vortices, 
there  are  three  dimensionless  parameters:   N  (discrete), 
(Ka)  (continuous)  and  {y  /y)    (continuous)  where  N  >^  2, 
0  <_  (Ka)  <  00  and  -  oo  <  {y  /y)    <  oo  .   In  this  paper  we 
only  consider  the  motion  in  an  unconfined  space  although 
flow  within  (or  outside  of)  a  finite  closed  boundary  or 


bounded  by  a  seml-inf Inlte  wall  would  be  of  Interest*; 
the  study  of  such  flows  are  harder.   We  also  restrict 
ourselves  in  §3  to  the  usual  linearized  stability  problem, 
i.e.  a  normal-mode  type  of  analysis  for  small  perturbations 
of  the  position  (x,  ,y,  )  of  the  vortices.   Because  of  the 
number  of  relevant  parameters  this  much,  combined  with  the 
non-linear  numerical  calculations,  yields  a  considerable 
amount  of  information.   Other  interesting  stability  pro- 
blems are  1)  perturbations  of  K  and  2)  perturbations  of 
y,;    these  problems  may  be  considered  in  the  future. 

We  make  a  few  remarks  about  the  possible  applicability 
of  our  results  to  the  description  of  actual  fluid  motions. 
Since  the  rectilinear  geostrophlc  vortices  are  derived  from 
a  tangent  plane  approximation  of  atmospheric  motion  over 
the  earth,  some  discretion  must  be  used  in  trying  to  apply 
them  directly  to  the  study  of  global  motions  in  which 
there  are  large  amplitude  meridional  (north-south) 
motions.   Although  the  original  motivation  for  these  studies 
has  come  from  meteorological  considerations,  we  believe  our 
results  may  be  relevant  in  the  description  of  other  fluid 


*For  example,  Onsager  [9]  has  studied  the  problem  of 
clustering  of  N  logarithmic  vortices  of  arbitrary  strength 
7.  moving  inside  of  a  closed  boundary. 


mechanical  phenomena*  In  which  the  forces  of  gravitation  and 
rotation  are  In  approximate  balance.  The  study  of  kinematic 
Instability  may  be  relevant  to  the  turbulence  problem. 

In  a  later  report  we  shall  present  the  results  of  our 
nonlinear  numerical  computations,  which  will  clarify  the 
range  of  validity  of  the  linearized  stability  analysis. 


^Recently  geostrophlc  vortices  have  been  used  in  the  study 
of  superconductivity  [10]. 


2.   Equations  of  Motion  and  Equilibrium. 

For  this  (N  +  1)  -  vortex  problem,  the  positions  of  the 
N  circle  vortices  are  best  described  In  polar  coordinates 
(r.,0. )  and  the  position  of  the  center  vortex.  In  cartesian 
coordinates  (x  ,y„)-   The  non-linear  (kinematic)  motion  of 
the  k   circle  vortex  In  polar  coordinates  (r,  ,0,  )  satisfies 


(2.1)    r^=   -  K' 


N 


1=1 


ti 


.   ,    ^l^Pkl)  1 

sm  0,  .  

^^    Pkl   i 


2  %(Pko) 

-  7  K^  X   sin  9.     -   y      cos  0,  )   ^-  ^" 
'o^o      k^'o      k 


'ko 


(2.2) 


N  r 


«.  =  ^-^  2 


k  k 


1^  L 
1=1 


(r,  -  r,  cos  0,  .  ) 
k    1      ki' 


Ki(p 


^ki  J 


+  7  K^  r,  -  (x   cos  0,  +  y   sin  0,  ) 
'oj^k    ^o      k   ■'o      k' 


n  ^i(pko) 


ko 


where  (')  means  time  derivative;  K^  Is  the  Bessel  function 
of  the  second  kind  (first  order);  and 


p,  .  =  K(rf  +  r?  -  2r,  r.  cos  0,  .  )  ^^^  ,   0,  .  =  ( 0,  -  0.  ) 
^ki    ^  k    1     k  1      ki'     '    ki    ^  k    i' 


^ko 


=  "^f-k  -  2^i.(x.  c°s  ^,.  +  y.  si>^  ^v)  +  (^?  +  y?) 


-11/2 


k'  o 


k 


k' 


The  non-linear  motion  of  the  center  vortex  In  cartesian 


coordinates  i^QtY^)    is  described  by 


(2.3) 


N   r 


.^=-K' 


(Jq  -  ^^   sin  e^) 


Kn  (p  .  )  -1 
Poi 


(2.4) 


o  N  r  K,  ( p  . )  ~i 

y  =  K^  r~  (x   -  r.  cos  0.)  — — 2^ 


where 


oi 


r .  -  2r.  (x  co£ 

1       1^  o 


e.    +  y      sin  0.  )  +  (x^  +  y^)  -^/^ 

1    "^  O        1        O    "^  o 


The  following  dimensional  normalizations  have  been  made: 
lengths  are  normalized  by  the  circle  radius  a;  time  Is 
normalized  by  (27ra  /y) ;    and  the  center  vortex  strength  7 
Is  normalized  by  7.   The  2(N  +  1)  Eqs .  (2.1)  -  (2.4)  state 
that  the  velocity  of  a  particular  vortex,  say  the  k   one. 
Is  the  sum  of  the  field  velocities  of  the  remaining  vortices 
evaluated  at  the  k   vortex  position,  each  field  velocity 
being  directed  normal  to  the  distance  vector  from  each 
vortex  to  the  k   vortex;  changing  the  sign  of  the  strength 
7   (or  7)  reverses  the  direction  of  rotation  of  the  vortex. 
In  this  dynamically  balanced  system,  each  vortex  moves 
slowly,  like  a  massless  particle,  in  the  velocity  field  of 
the  other  vortices. 


2.1   Equilibrium  Motion. 

For  N  vortices  of  the  same  strength  7  equally  spaced 
on  a  circle  and  one  vortex  of  strength  7  at  the  center  of 
the  circle,  the  circle  vortices  move  with  uniform  angular 
velocity  O  around  the  center  vortex.   For  this  equilibrium 


7 


configuration  the  circle  vortices  satisfy  the  following 
conditions : 


r^^  =  1  ;    r^^  =  0 


0^  =  w^  ,    ^k  "  "k  "  ^  '   const, 


where 


Wj^  =  Ofc  +  w^  , 


and  the  superscripts  denote  the  initial  equilibrium  values 


w 


.0  _  2vk 
k    N 


The  center  vortex  satisfies 


X   =y   =0:   X   =y   =0 
o   -^o     '    o   "^  o 


Inserting  these  conditions  in  Eqs.  (2.1)  -  (2.4)  we  obtain 
for  the  circle  vortices 


(2.5) 

and 
(2.6) 


N   r 


1-1 


sm  w 


■^i(\i'  -I 


ki  a 


kl 


=  0 


"  =  %  +  ^o"^o 


O.  =^I 


N 


1^  L 
1=1 


kl  1  ^  kl ' 


0)^  =  KK-^(K) 
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where  w,  .  =  (w,  -  w.)  =  (w,  -  w.  )  and  a     .    = 


2(1  -  cos  w,  .  ) 
.  ki 


1/2 


The  center  vortex  equations  yield 


the  conditions 


(2.7) 


N  N 

5    sin  w.  =  0  ,    \        cos  w.  =  0 

fci         ^  to         ^ 


For  example,  the  equilibrium  configuration  for  N  =  8  is 
shown  in  Fig.  1. 


3-   Linearized  Stability  Analysis  of  N  Vortices  on  a 
Circle  and  One  Vortex  at  the  Center. 

The  linearized  equations  of  motion  are  obtained  in 
the  usual  way.   We  make  a  perturbation  expansion  of  the 
position  vector  (r.,e.;  x  ,y  )  on  the  equilibrium 
solution;  the  perturbation  parameter  is  e  <<  1.   For  the 
k   circle  vortex 

(3.1)    ^k  =  ^  +  ^4^^  "^  °(^^^ '  ^k  =  ^^k"^^  +  °(^^) 


where  the  equilibrium  angular  velocity  O  =  w,  is 
given  by  (2.6).   For  the  center  vortex 

(3.3)  x^  =  ex^^^  +  0(e2)  ,   x^  =  ei^^^  +  0{e'^) 

(3.4)  y^  =  ey^^^  +  0(6^)  ,   y^  =  ey^^^^  +  0(8^)  . 


The  expanded  solution  (3-1)  -  (3-4)  is  put  into  (2.1)  -  {2A) , 

retaining  only  first  order  terms  in  e.   The  presence  of  the 

center  vortex  introduces  periodic  coefficients  into  the 

resulting  system  of  linear  differential  equations  for 

(r''  ^^l■      'f    y^ry      ' ^Q      ''      However,  this  apparent  difficulty 

can  be  eliminated  by  transforming  the  equations  to  a  frame 

of  reference  rotating  with  the  equilibrium  angular  velocity  O, 

x^^^  =  x|^^^  cos  Ot  +  y^^^  sin  Ot 

y^-^^  =  -  x^^^  sin  Ot  +  yj^^^  cos  -nt  . 
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The  resulting   system  of  2(N  +  1)    linear  differential   equa- 
tions  for   (r^     >^]^   '    ^0^    '^Q^    '    with  constant   coefficients 


are 


i-1 


i=l 


7.^.(xi^^sin  w°   -  ^i}^    cos  w°)    =  0 


O     C      Q 


ic         "  O 


n(l)      ^     fr..n         ,      „.     .     ^  Jl)      .    X^    c      .r!l)      4-    T^    A      .«(!) 


(3.6)        e;^^'   +  (O  +  C^  +  7,Djr^^^    4- 


i^k      ^i    i         ■   fe  ^'^i^i 


i=l 


i-1 


-  7.Dj3?i^)    cos   wj  +  5^^^)    sin  w°)    =  0 


o   o'    o 


^(1) 


N 


,(1)    ...      o 


(3.7)        x^-"^    +  D^  J-     r-'    sin  wT    -  cd 

O  O    'T_.,         1  10 


N 


5^^ 


(1)  o 

;         cos   w. 

1  1 


-[' 


N 


(o)     +  D    )    >        sin  w?   cos   w? 
-  i-1  ^  i 


~(1) 
^o 


N 


2    o 


(o)^   +  D^)   \        sin     wr    -  NcD     + 


O 


e'  =  ° 


N 


(3.8)        ¥^^    -  D^  :^  rjl)    cos  w°       o.^  TZ  ^1^^    sin  w? 
o  o  ■^■^T-      1  10  'f~-=r      i  1 


(oo     +  D^)   "S       cos   ^w?   -    Noo     +  O^x^^ 
_      O  O'    /^~-  1  o  Jo 


T-  N  -I 

+      (oc     +  D^)   5        sin  w?   cos   w° 

J_        O  O       <  _;'-|-  1  1 


e'  =  0 
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where 


-K' 


ki 


sin  w,  .  K  (cr,  .  ) 
ki   o^  ki 


ki    2 


^ki  -  "5- 


o 


(1  +  cos  w."".  )  K  (a,  .  )  + 
^         ki    o^  ki 


(1  -  cos  w,  .  )  K  (a,  .  )  +  ^ 

^  ki'   o^  ki'      ? 


^ki   J 


ki   -J 


,2  N   r 


^o  --- 


37k 
1=1 


ki 


(1  -  cos  W°.  )  K^(cT^i)  -  2  COS  w°. 


v^kiin 


2(1  -  cos  w 


kl' 


1/2     o    2ir    f,  .  . 

'    ^ki  =  -r  (^  -  ^) 


CD^  =  KK^(K)  ,    Dq  =  "^^ 


K  (K)  + 

o^  ' 


Ki(K)  1 


X" 


Eqs.  (3-5)  -  (3.8)  are  valid  for  N  >  2.   For  N  >  3  the 

center  vortex  equations  (3.7)  and  (3.8)  simplify  because  of 

Increased  geometrical  symmetry  expressed  by  the  additional 

conditions* 

N  N 

(3.9)   >   sin  2  w.  =  2  >   sin  w.  cos  w.  =  0  , 
1=1        ^1=1      ^      ^ 


N  N      2        ^ 

>   cos  2w.  =  ">   (cos  w.  -  sin  w.  )  =  0 

1=1       ^1=1       ^        ^ 


.  2 


*A.  Troesch  has  Informed  us  that  the  conditions  (2.7)  and 
(3.9)  SI'S  included  in  a  function  theoretic  result  due  to  Veska, 


12 


and  the  Identity 

N      p        p 
(3.10)  J        (sln^'w.  +  cos  w.)  -  N 


5.1  Computation  of  Eigenvalues. 

Since  the  linearized  equations  of  motion  (3.5)  -  (3.8) 
have  constant  coefficients  the  usual  exponential  stability 
analysis  is  possible.   The  displacement  (rJ'   >^i   ,  x  ,y  ) 
from  the  equilibrium  position  can  be  written  as  a  sum  of 
eigensolutions  in  the  form  exp.  (A^t)  where  Z  =   1,2, . . .2(N+1) 
The  motion  is  unstable  with  respect  to  the  equilibrium 
position  if  any  eigenvalue  Ao  has  a  positive  real  part. 

If  the  center  vortex  is  absent  (7  =0),  there  is  a 
high  degree  of  symmetry  and  the  eigenvalue  calculation  can 
be  reduced  to  solving  a  quadratic  equation.   For  7   =0, 
the  reduced  Eqs.  (3.5)  and  (3.6)  take  the  matrix  form 

A    B 


(3.11) 


=  ^i 


where  j  =  1,2,...N,  i  =   1,2,...2N  and  the  sub-matrices 
A,  B  and  C  are  each  N  x  N  right-circulant  matrices.   For 
example,  the  elements  of  A  for  N  =  4  are 
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A  = 


^11  ^12  ^13  ^1^ 

^14  ^11  ^12  ^13 

^15  ^14  ^11  ®12 

'l2  ^13  ^14  ^11 


/ 


The  eigenvectors  of  a  circulant  matrix  are  independent  of 
the  matrix  elements  and  are  given  by 


(1,  J,    a.2J,  ...,  co(^-^)J) 


where 


(D  =  e 


Then  the  sub-matrices  A,    B  and  C  all  have  the  same  eigen- 
vectors and  hence  commute;  their  respective  eigenvalues  are 

(3.12)  a^.  -  a^^  +  a^^ui^   +  ...  +  Sufo^^"'''  "^",  J  =  1,2,...,N 

(3.13)  bj  =  b-j_^  +  h^^a^^  +   ...    +  \^[^~'^'^^ 


(3.14)   Cj  =  c^-^  +  c-^2^^  + 


+  c 


llf 


(N-l)j 


The  eigenvalues  A.  in  (3. 11)  are  given  by 

1/2 


(3.15) 


A,  =  a  .  i  (b.c  .)• 


which  agrees  with  Stewart  [1]. 

In  the  general  case  when  the  center  vortex  is  present 


14 


(7o  ^   0),  the  clrculant  property  is  destroyed  and  the 
eigenvalue  calculation  for  the  resulting  2(N  +  1)  x  2(N  +  1) 
matrix  Is  carried  out  by  an  Iteration  procedure  devised  by 
B.  Parlett  [11].   The  [y ^  =   0)  case  has  been  used  to  test 
the  accuracy  of  this  Iteration  method. 


5.2   N  Vortices  on  a  Circle  (7   =0). 

This  case  was  considered  by  Thomson  [9]  for  logarthmlc 
vortices  (K  =  0)  and  Stewart  [1,2]  for  the  geostrophlc 
K^-vortlces  (K  >  0). 

A   physically  significant  transition  region,  depending 
on  both  K  and  N,  Is  found  separating  stable  and  unstable 
configurations.   In  Fig.  2  the  variation  of  the  maximum 
real  eigenvalues  with  respect  to  K  for  2  <  N  <  9  is  shown. 
For  N  =  6  the  circle  configuration  is  stable  for  0  <  K  <  1.288 
and  unstable*  for  K  >  1.289.   For  N  >  7  the  circle  configuration 
is  unstable  for  all  K;  except  K  =  0  is  a  neutral**  stability 
point  for  N  =  7.   For  2  ^  N  <_  5  the  circle  configuration  is 
stable  for  all  K. 

In  Fig.  3  the  variation  of  the  Imaginary  eigenvalues 
(normal  or  elgenfrequencies )  with  respect  to  K  for  N  =  6  is 


*Thls  corrects  Stewart's  result  that  N  =  6  is  stable  for 
all  K. 

**Thls  corrects  Thomson's  slightly  positive  eigenvalue. 
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shown.   The  equilibrium  frequency*  O  is  marked.   All  of  the 
normal  frequencies  decrease  with  increasing  K.   However, 
it  is  physically  more  meaningful  to  compare  all  frequencies 
with  O  (or  O).   Since  in  general  O  decreases  more  rapidly 
than  the  other  eigenfrequencies  with  increasing  K,  vortices 
with  larger  K  behave  like  particles  with  harder  potentials 
in  dynamical  systems.   In  addition  there  are  double  zero 
eigenvalues  for  all  values  of  N  and  K.;    these  are  discussed 
in  §3-6. 

5.3   (N  +  1)  Vortices  {y ^  =   1). 

This  case,  with  N  vortices  on  a  circle  and  a  center 
vortex  of  the  same  strength  as  each  of  the  circle  vortices. 
Illustrates  some  of  the  main  features  of  the  general  case 
(7  _>  0).   Comparison  with  the  case  (7   =  0)  is  particularly 
instructive. 

Similar  to  the  case  without  a  center  vortex,  a  physically 
significant  transition  region  is  found  separating  stable  and 
unstable  configurations  for  7  =1.   In  Fig.  4  the  variation 
of  the  maximum  real  eigenvalues  with  respect  to  K  for 
2  <  N  <  10  is  shown.   For  N  =  8  the  (N  +  l)-vortex  configuration 
is  stable  for  0  <  K  <^  1-597  and  unstable  for  K  >  I.598.   For 
N  >  9  the  configuration  is  unstable  for  all  K;  except  K  :=  0 


*Thomson  calculated  an  extraneous  eigenfrequency  larger  than 
Q  for  K  =  0  and  N  =  5  and  6  due  to  misusing  the  second  of  the 
t^o  linear  integral  invariants  (3.17)  and  (3.I9). 
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is  a  neutral  stability  point  for  N  =  9.   Comparing  with  Fig.  2 

we  see  that  N  =  8,  9>  10  for  y      =1  are  comparable  to 

N  =  6,  7,  8  respectively,  for  7   =  0.   That  is,  the  presence 

of  a  center  vortex  (7   =  1)  allows  the  presence  of  two 

additional  circle  vortices  without  appreciably  changing  the 

stability  characteristics.   However,  for  N  =  2,  this 

(N  +  1) -configuration  is  unstable  for  all  K;  for  N  =  3,  this 

configuration  is  unstable  for  0  jf,  1^  _f,  3-^5  and  stable  for 

K  >  3«50'   Thus  for  an  intermediate  number  of  circle  vortices 

this  (N  +  1) -configuration  is  stable  for  all  K,  i.e.,  for 

N  =  ■4-,  5,  6  and  7. 

Comparison  with  7   =0  implies  that  the  instability 
for  large  N  (>  6)  means  an  unstable  circle  vortex.   The 
instability  for  small  N  (<  4)  means  an  unstable  center  vortex; 
this  is  clearly  seen  by  geometrical  considerations  for  N  =  2 
and  7  — >  0  (cf.  Fig.  6,  K  =  0).   These  characteristics 
become  more  apparent  in  the  following  sub-sections. 

In  Fig.  5  the  variation  of  the  normal  frequencies 
with  respect  to  K  for  N  =  8  is  shown.  The  equilibrium 
frequency  O  =  Q  +  J J^cSy ^   =  1)  Is  marked. 

3.4   (N  +  1)  Vortices,  General  Case  (-00  <  7   <  00 ) . 

We  are  now  in  a  position  to  describe  the  general  case 
where  N_>2,  0_<K<oo  and  -  00  <  7   <  00  .   The  simplest 
approach  is,  first,  to  study  the  special  case  of  logarithmic 
vortices  (K  =  0)  which  is  representative  of  the  7   variation 
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of  the  eigenvalues  for  all  K. 

In  Fig.  6  the  variation  of  the  maximum  real  eigenvalues 

with  respect  to  7  for  K  =  0  and  2  <  N  <  10  is  shown.   Two 
^       '  o  —   — 

distinct  families  of  real  eigenvalues  are  evident:   one 

family,  which  yields  (for  different  N)  stable  motions  for 

7   greater  than  the  neutral  stability  point,  is  pertinent 

for  the  N  circle  vortices;  and  the  other  family,  which 

yields  (for  different  N)  stable  motions  for  7   less  than 

the  neutral  stability  point,  is  pertinent  for  the  center 

vortex.      Note  that  for  N  =  2  the  unstable  range  includes 

7   =0,  i.e.  unstable  center  vortex.   For  example,  for  N  =  4 

the  region  of  stable  motions  is  -  O.5  <  7   <  2.25;  fo^ 

7   <  -  0.5  a  circle  vortex  is  unstable  and  for  7  >   2.25 
'o—  'o— 

the  center  vortex  Is  unstable.   In  Fig.  7  the  variation  of 
the  normal  frequencies  with  respect  to  7   for  K  =  0  and 
2  <_  N  <_  6  is  shown.   The  equilibrium  frequency  O  varies 
linearly  with  7  as  well  as  with  N. 

The  linearized  stability  of  (N  +  1)  vortices  over  the 
entire  range  of  the  three  parameters  N,  K  and  7  are 
summarized  in  two  figures,  Fig.  8  and  Fig.  9=  1)  Fig-  8  shows 
the  neutral  stability  curves  for  the  circle  vortices  where  the 
stable  region  is  to  the  right  of  each  curve  and  2)  Fig.  9 
shows  the  neutral  stability  curves  for  the  center  vortex 
where  the  stable  region  is  to  the  left  of  each  curve.   Both 
figures  must  be  superposed  in  order  to  determine  the  regions  of 
stability  (or  instability)  for  each  N.   For  N  =  2  and  J>   there 


18 


is  a  seml-inflnlte  range  of  stability  for  large  negative 

7  ,  I.e.  the  circle  vortices  are  stable  for  all  7  .   For 
o  o 

N  _>  4  the  vortex  configuration  is  stable  in  a  finite  range 
of  7  and  the  size  of  this  range  increases  rapidly  for 
large  N;  this  is  shown  in  Table  I.  for  K  =  0. 


3.5   (N  +  1)  Vortices  with  a  Fixed  Center  Vortex. 

The  vortex  configurations  studied  in  the  previous 
sections  were  free  vortices  without  external  constraints. 
In  this  section  we  consider  the  effect  of  holding  the 
center  vortex  fixed  at  the  center  instead  of  being  free  to 
move.   Then  Eqs.  (5.5)  -  (3«8)  are  modified  by  the  addi- 
tional conditions. 

^(1)   -^(1)    n  .  ^  ~(1)    '^(1)    n 
x^    =  y    =0  and  x^  '  =  y^  '  =  0  . 
Q     -^G  o     •'o 

The  number  of  degrees  of  freedom  is  2N  instead  of  2(N  +  1) 
for  the  free  case.   Eqs.  (^-T)  and  (3-8)  are  identically 
zero;  and  Eqs.  (5.5)  and  (3.6)  now  have  the  circulant 
property  discussed  in  §3.1  for  the  case  (7   =  0)  without  a 
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Table  I 
Range  of  Exponential  Stability  for  Free  Center  Vortex  Case,  K  =  0 


N 

Lower  Stability  Limit 

Upper  Stability  Limit 

2 

7o  =  -00 

7o  =  -1.25 

5 

-00 

1 

4 

-.5 

2.25 

5 

-.5 

4 

6 

-.25 

6.25 

7 

0 

9 

8 

•  5 

12.25 

9 

1 

l6 

10 

1.75 

20.25 

11 

2.5 

25 

12 

3.5 

30.25 

15 

4.5 

36 

14 

5-75 

42.25 

15 

7 

49 

A  circle  vortex  is  unstable  for  y      less  than  the  lower 
stability  limit. 

The  center  vortex  is  unstable  for  y     greater  than  the 
upper  stability  limit. 
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center  vortex.   Again  the  eigenvalue  calculation  can  be 
reduced  to  solving  a  quadratic  equation,  the  only  differ- 
ence from  the  case  (7   =  0)  being  a  modification  of  right 
diagonal  terms  In  the  submatrlx  C,  I.e.  the  coefficient  of 
r^-*-'  Is  (O  +  Cq  +  7q^q)    Instead  of  (O^  +  C^)  for  7=0 
(cf.  [2]). 

The  main  results  are:   1)  for  N  >^  4  the  Instability 
characteristics  of  the  circle  vortices  are  the  same  for 
fixed  and  free  center  vortex.  I.e.  same  real  eigenvalues, 
2)  for  N  =  2  and  3  the  circle  vortices  have  a  semi-Infinite 
unstable  range  for  negative  7  :  for  example  for  K  =  0  the 
circle  vortices  are  stable  for  N  =  2  over  the  range 
7   >  -  0.25  and  for  N  =  3  over  the  range  7   >  -  O.5 
(cf.  Table  I  for  free  center  vortex)  and  3)  the  equilibrium 
frequency  O  is  not  an  elgenfrequency  if  the  center  vortex 
is  fixed;  and  2(N  -  2)  of  the  eigenvalues  are  the  same  as 
the  eigenvalues  for  the  free  (center  vortex)  case,  including 
the  double  zero  roots.   This  is  an  example  in  which 
constraints  have  introduced  some  instability. 
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3.6  Multiple  Elgenfrequencies 

Over  the  entire  range  of  the  parameters  N,  K  and  7 
the  zero  eigenvalue  is  at  least  a  double  root;  on  the 
neutral  stability  curves  shown  in  Figs.  8  and  9  there  are 
quadruple  zero  roots.   In  addition  there  are  a  countable 
infinity  of  multiple  non-zero  eigenfrequencies,  especially 
for  K  =  0.   By  standard  methods  [12]  we  can  demonstrate 
that  the  multiple  non-zero  eigenvalues  do  not  produce 
solutions  of  the  form  t^~   •  e   ,  where  |j,  is  the  multiplicity. 
However  the  multiple  zero  roots  require  some  attention. 

There  are  two  linear  invariants  associated  with  the 
system  of  equations  (3.5)  to  (3.8)  inclusive.  Either  by 
summing  the  N  equations  {3- 5)  or  by  linearizing  the  non- 
linear integral  invariant  (A5.2),  we  get 

^   .(1)  _ 


(3.16)         i_  K     =  0 

k-l  ^ 
and 

N   (-1) 
(3.17)  y        r;  '  =  const.  =  0 

k=l   ^ 

The  constant  in  (3.1?)  can  be  set  equal  to  zero  by  an 
appropriate  choice  of  the  initial  equilibrium  radius. 
Either  by  summing  the  N  equations  (3-6)  or  by  linearizing 
the  non-linear  integral  invariant  (a6.1),  we  get 
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N   .,-,,      ,  N 


(3.18)  ytel^^    --(t-rl^^)[o.C      .,D   ./-C  1 

k=l  ^      ^fa  ^   ^  L     o   ^o  o  4^   ^kl_ 


i=l 
const.  =  0 


The  constant  in  (3.l8)  is  zero  only  If  the  constant  in 
(3.17)  is  zero.   Then  (3.I8)  implies 


(3.19)  XI  ^k  "  ^°^s^- 


Since  the  second  linear  invariant  (3.I9)  is  a  consequence 
of  the  first  linear  invariant  (3-17),  the  constant  in 
(3.19)  cannot  be  set  equal  to  zero  for  all  possible  initial 
conditions. 

For  7o  =  0  the  linear  invariants  (3.I7)  and  (3.I9) 
can  be  associated  with  the  double  zero  eigenvalues, 
eliminating  unstable  solutions  of  the  form  t^""*"  where  p.  is 
the  multiplicity.   But  at  a  neutral  stability  point  such 
as  K  =  0  for  N  =  7,  the  quadruple  zero  eigenvalues  yields 
a  linearly  time-increasing  solution. 

The  linear  invariants  (3.17)  and  (3-19)  are  apparently 
valid  for  7^  7^  0;  but  since  these  invariants  are  independent 
of  the  position  of  the  center  vortex,  we  suspect  that  a 
linearized  description  of  the  motion  may  not  be  strictly 
valid.   If  these  invariants  cannot  be  invoked,  the  double 
zero  roots  yield  an  algebraic  instability  proportional  to 
t  over  the  entire  range  of  parameters  (except  for  7   =0). 
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The  implication  is  that  the  center  vortex  motion  is 
Inherently  non-linear.   Thus  we  must  study  the  non-linear 
initial  value  problem  in  order  to  determine  the  range  of 
validity  of  the  linearized  stability  analysis.   In  deriving 
the  linearized  equations  (3«5)  to  (3.8),  some  peculiarities 
of  the  center  vortex  were  already  evident:   first,  mixed 
cartesian  and  polar  coordinates  were  needed  to  describe  the 
motion  of  the  center  and  circle  vortices;  and  second,  the 
center  vortex  terms  introduced  periodic  coefficients  which 
fortuitously  could  be  made  constant  by  transformation  to 
an  appropriate  rotating  coordinate  system. 
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Appendix  I.   Integrals  of  the  Motion. 

The  number  of  Integral  Invariants  or  constants  of  the 
motion  Is  the  same  as  for  the  equations  of  motion  of 
logarithmic  vortices  (cf.  Lamb  [I3]).   Likewise  the 
kinematic  equations  of  motion  for  geostrophlc  vortices  can 
be  written  In  a  more  symmetric  form  than  those  given  In 
§2 

''k     k 
where  the  Klrchhoff  fixnctlon  W  for  geostrophlc  vortices  Is 

1   N 
(A2)  W  =  ^  y~  y.y  .    K  (p.  .) 

47r  4^  ''i^j   o^^ij' 

and  p.  .  =  K[(x.  -  X  )^  +  (y   -  y.)^]""-/^;  for  K  ==  0,  W  reduces 
to  the  original  logarithmic  expression  of  Klrchhoff.   In 
polar  coordinates  (Al)  becomes 


The  Integral  Invariants  are: 

1)   W  =  0  and  W  Is  an  Integral  Invariant.   This  follows  from 
the  symmetry  of  W  In  1  and  J;  and  because  W  Is  not  an 
explicit  function  of  time  and  Is  only  a  function  of  the 
distance  between  vortices  p.  .. 
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2)   From  (Al)  and  the  translatlonal  invarlance  of  W 


(M.2)         IZr^yk  =  -  IZ  If-  =  ° 

k^  ^  ^     k=l  "^^k 


and,  for  example,  the  centrold  is  an  integral  invariant 

N                          N 
(A4.3)   jj =  const.,       Tj =  const. 

N  N 

if  >   7,  ?^  0;  for  >   7,  =  0,  the  numerators  are  still 
k=l  ^         k-1  ^ 

constants  but  the  physical  interpretation  is  less  apparent. 


3)   From  (A5.I)  and  the  rotational  invarlance  of  W 

(^5.1)  II  Vk^k  =  21:  3^=  0 

k=T  ^  ^  ^   k=T  ^^k 

and 


N 
(A5.2)  3   7,  rf  =  const. 


The  invariant  (A5.2)  for  this  kinematic  system  is  comparable 
to  the  momentum  invariant  in  dynamical  systems. 
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4)   From  (A3. 2) 


r, 


k=l  '^  '^  '^     k=l  ^  ^^k 


N 


l,k=l 


=  const. 

since  the  right  hand  side  of  (a6.1)  has  the  same  syTimietry 
properties  as  W.   For  K  =  o,  (A6.1)  reduces  to  the  classical  result, 

i,k=l 
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Fig.  1  Equilibrium  Configuration  for  N  =  8 
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vortex  of  strength  7  at  the  center  of  the  circle.   One  length  scale 
K~   is  associated  with  each  of  the  vortices,  the  limiting  case  K.   =  0 
being  the  classical  logarithmic  vortex.   The  exponential  stability 
of  this  equilibrium  configuration  is  studied  over  the  range  of  the 
three  parameters:   N  ^  2,  0  ^   (Ka)  <   co  and  -od  <  iy^/y)    ^   00.   The 
need  for  a  supplementary  nonlinear  study  is  discussed. 
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